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(Dated: May 31, 2010) 

In this comment we show that there is a direct connection between coherent exchange of energy 
among light and matter and the emission spectrum of a microcavity quantum dot system as modeled 
in Phys. Rev. B 79, 235325 (2009) by F. P. Laussy, E. del Valle, and C. Tejedor. To do so, we 
show that in their model the necessary and sufficient conditions for having eigenvalues with non-zero 
imaginary parts in the propagator of the bare mode populations, are the same as for having strong 
coupling in the emission spectrum. This amounts to saying that, whenever there is strong coupling 
there will be oscillating frequencies in the dynamics of the populations. These conditions are valid 
both for the case where matter is treated as bosonic or fermionic, in the spontaneous emission case. 

PACS numbers: 42.50.Ct, 78.67.Hc, 42.55.Sa, 32.70.Jz 
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Laussy and coworkers ^"'^ studied the effects of consid- 
ering incoherent pumping of both excitons and photons 
in a microcavity quantum dot system. In particular, they 
clearly show the nontrivial effects of such pumping mech- 
anisms in the emission spectrum of the system. The 
system is modeled using the following quantum master 
equation (in units where h — \): 
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where H = J2c=a b ^cC^'^ + 9 ("^^^ + '^^^) • '^^'^^ ^-^^ pho- 
tonic boson operators and 6, 6^ can be either bosonic or 
fermionic matter operators. The quantities Wc are the 
bare energies of the light and matter, g is the coupling 
constant between them, 7t, are the decaying rates of the 
cavity and the emitter and Pc are the rates at which they 
are being incoherently pumped. 

For the case where the exciton is modeled as a boson 
Laussy et. al. define a criteria for having strong coupling 
(SC)i: 
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where 4r± = r„ ± Fb = (7, - P^) ± (7^ - Pb). We 
will show, that the same criterion of the above equation, 
is a necessary and sufficient condition for having oscil- 
latory frequencies in the propagator of the bare mode 
populations. For convenience we re-write equation (12) 
of reference [l| as follows: 



(3) 




where Ua — {a'^a), ni, = (&^&), « = ^{{a''b)), (3 = 
3((a^6)), and A = Wq — ojf,. The above equation can be 
written in a more compact form as: ^w(t) = f + Aw(i), 
where A corresponds to the coefficient matrix, w(i) is 
the vector which contains the single time mean values 
of interest, i.e., w(i) = {na{t),nb{t), (3{t),a{t))'^ and 
f = {Pa, Pb, 0, 0)-^. The formal solution of this equation 



w(i) = e^*w(0) + A-i (e^* - l) f. 
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where e^* is the propagator of the last equation. The 
eigenvalues of the matrix A are given by: 
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In resonance condition (A = 0) they simplify to: 

A(A) = {-2F+, -2r+, -2r+ + 2ii?o, -2^+-2^i?o}, (6) 



where Rq = J g"^ — F_ . Notice that if i?o is a positive 

number, the condition ([2|) is automatically fulfilled, i.e., 
light and matter are in the SC regime. Then, if i?o > 0, 
there will be imaginary frequencies in the propagator that 
will lead to oscillations in the populations. As a conse- 
quence of choosing as variables the real and imaginary 
parts of the coherence (a^6), the propagator in resonance 
condition is given by a block diagonal matrix, with form: 
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Notice that the variable a decouples from equation ^ 



and undergoes exponential decay a{t) = 
The square matrix M3x3(i) is given by: 
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The above equation explicitly shows that in SC the popu- 
lations undergo oscillations with frequencies proportional 
to 2i?o. 

In this paragraph, we re-examine the results in section V 
of ref. PI . In figure 6 of the same article they present the 
dynamics of the average photon number n^ for different 
initial conditions. In particular, they claim that in strong 
coupling there is no oscillation in such observable for the 
initial condition of one upper polariton with the param- 
eters A — Pa — Pb = 0, ^a = 1-95 and jb — O.lg. When 
one plugs their parameters into the propagator, equation 
([5]), one sees that there are oscillating frequencies, yet 
the dynamics seems to be simply decaying in the popu- 
lations Ua and Hb- When their difference is plotted it is 
seen that it actually oscillates showing that indeed both 
terms contain imaginary frequencies (see the left panel of 
fig. [T]). One might wonder why the oscillations are not 
clearly visible in the populations. The reason is that the 
spectrum of A contains two purely real eigenvalues that 
will hinder the other two which are complex. This is fur- 
ther clarified when one takes the Fourier transform of the 
propagator and applies it to a given initial condition. In 
the right panel of fig. [T]we plot the Fourier components of 
Ua for the parameters of the left panel of the same figure. 
It is clearly seen how the widened Lorentzian around zero 
hides the rest of the contributions at ±2ii!o ~ ±1.786065. 
The reason why Ub — Ua clearly shows oscillations, with 
period T = 2tt/{2Ro), is precisely that the zero frequency 
component cancels out after the substraction. 
Summarizing the above discussion, we have shown that 
for the quantum master equation ([1} the following state- 
ments are equivalent: 

• The propagator of the bare mode populations has 
imaginary frequencies. 

• The system, as modeled by equation ([T|), is in SC 
regime. 

This implies that, whenever oscillations are observed in 
the dynamics of the bare mode populations the system 
will be in SC regime. In the case where one models the ex- 
citon with a fermionic operator [b — a), the spectrum of 
photoluminescence cannot be obtained analytically when 
the effects of incoherent pumping are present. When 
Pa = P(, = one can obtain exactly the eigenvalues of 
the dynamical matrix involved in the equations of the 
quantum regression theorem, del Valle et. al. define nth 
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FIG. 1: (left panel) Time evolution of the populations of 
the photons Ua, excitons rib and their difference for the pa- 
rameters A = Pa = Pb = 0, 7a = 1.9(7 a-nd 71, = O.lg. The 
initial condition was an upper polariton, {a^6)(0) = na(0) = 
711,(0) = 1/2. (right panel) The dashed lines are the compo- 
nents of the Fourier transform of the signal na{t). The solid 
black line is the sum of the components of the signal, i.e., the 
Fourier transform of Ua (i) . 



order SC as the following condition^.: 

g>\r-\/V^ 



(9) 



Here, we first show that the condition above is also 
necessary and sufficient to have imaginary frequen- 
cies, in the dynamics of the populations and co- 
herences in the same model. To do so, we use 
the notation defined in [J|. We arrange the pop- 
ulations and coherences of the density matrix p 
as follows: u == {pgo.,go, PG1..G1, Pxo.xo, l3i,ai, . . . , 

PG7i,Gn,PX7i-l,Xn-l,Pn,Ctn, ■ ■ ■) , Oin = ^{pGn,X7i-l) , 

Pn = 3(pGn,Xn-i) and Pim.jn = {i,m\p\j ,n) . Where 
i, j are either G (ground state) or X (exciton state) and 
n, m are integers representing the number of photons in 
the Fock state. 

The equation of motion of u is 



— u = Bu. 

dt 



(10) 



The structure of the matrix B will consist of block diag- 
onal terms of sizes 1x1 (only one block, corresponding 
to the vacuum pgo,go)j 4x4 and off diagonal terms over 
the diagonal of the matrijsi. This implies that to reduce 



the matrix B to an upper triangular matrix it is neces- 
sary to rotate each diagonal block and the blocks above 
it. Once the matrix has an upper triangular form their 
eigenvalues are simply given by the elements of the diag- 



onal. Summarizing, the eigenvalues of the whole matrix 
B are simply the eigenvalues of the blocks 1x1 and 4x4. 
The structure of the blocks 4x4 corresponding to the 
dynamical equations of the nth excitation manifold is: 
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The eigenvalues of the above matrix are given by; 
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Notice that all equations derived for fermionic matter 
reduce to the boson matter model when n = 1 and 
Pa = Pf, = 0. However this correspondence breaks down 
with finite pumping or for higher excitation manifolds. 



In order to have c ± \d 



A=o 



< it is necessary that 



g > \ia — 7f)|/(4-\/n) — iT-l/v^- That is, the condition 
for having nth order strong coupling is the same as for 
having imaginary eigenvalues in the propagator of equa- 
tion PII)) . 

We would like to emphasize that the emission spec- 
trum of the system is a robust quantity, whereas the 



dynamics of the populations, in general, depends on the 
initial conditions. Yet, the propagators of the dynamical 
systems ([3]) and (|10l) are also robust features that show a 
similar dependence on the parameters that the emission 
spectrum presents. 

It is important to be able to correlate the dynamical 
regimes (SC or weak coupling, WC) with the dynamics 
of the populations and coherences. For instance, the 
study of the entanglement between the exciton and 
photonic subsystems will require the knowledge of the 
elements of the density matrix. In the SC regime one 
will expect entanglement sudden death and re-birth 
while in WC one will see how the entanglement either 
vanishes in a finite time or approaches asymptotically 
to zero^. Finally, we would like to emphasize that the 
connections found between population dynamics and 
the emission spectrum of the system are particular to 
the model given by equation ([T}. In general, one cannot 
expect them to be related, and one can find models 
where the frequencies of the propagator and the first 
order correlation function depend in a different way on 
the system parameters^"—. 



Electronic address: nquesada@pegasus. udea.e du. co| 
Electronic address: pcardenas@fisica.udea. edu.co| 
Electronic address: banghelo@fisica.udea. edu.co| 
F. P. Laussy, E. del Valle, and C. Tejedor, Phys. Rev. B 79, 

235325 (2009). 

E. del Valle, F. P. Laussy, and C. Tejedor, Phys. Rev. B 79, 

235326 (2009). 

F. P. Laussy, E. del Valle, and C. Tejedor, Phys. Rev. Lett. 
101, 083601 (2008). 

J. I. Perea, D. Porras, and C. Tejedor, Phys. Rev. B 70, 



115304 (2004). 

C. A. Vera, N. Quesada M, H. Vinck-Posada, and B. A. 

Rodriguez, J. Phys.: Condens. Matter 21, 395603 (2009). 

E. del Valle, Phys. Rev. A 81, 053811 (2010). 

N. Quesada, H. Vinck-Posada, and B. A. Rodriguez, 

arXiv: 1005.2318 (2010). 

C. A. Vera, H. Vinck-Posada, and A. Gonzalez, Phys. Rev. 

B 80, 125302 (2009). 



